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In this letter we review the separate universe approach for cosmological perturbations and point out 
that it is essentially the lowest order approximation to a gradient expansion. Using this approach, 
one can study the nonlinear evolution of inhomogeneous spacetimes and find the conditions under 
which the long wavelength curvature perturbation can vary with time. When there is one degree 
of freedom or a well-defined equation of state the nonlinear long wavelength curvature perturbation 
remains constant. With more degrees of freedom it can vary and this variation is determined by 
the non-adiabatic pressure perturbation, exactly as in linear theory. We identify combinations of 
spatial vectors characterizing the curvature perturbation which are invariant under a change of time 
hyp ersurf aces. 



Until quite recently, a standard tool used in inflationary calculations was the conservation of the curvature pertur- 
bation £ on superhorizon scales. This made inflationary predictions insensitive to what was going on between horizon 
crossing and horizon reentry and thus could lead to very robust predictions. However, it was realised Q that in the 
presence of more degrees of freedom £ can vary on superhorizon scales due to the presence of a non-adiabatic pressure 
perturbation. This fact was shown in a very simple manner in 1] using only energy momentum conservation. The 
rule for a varying or non-varying £ was derived using linear perturbation theory. The extent to which the rule is 
true when nonlinearities are taken into account remained largely unanswered, although one would expect intuitively 
that it would still hold. As far as we are aware, there are three references where nonlinear conservation of a variable 
connected to the linear £ has been mentioned 0-0. Q| and, in all cases, their results refer to single field inflationary 



Dealing with nonlinearities in general relativity is in general a very difficult issue. Yet there exists an approxima- 
tion which makes the problem tractable and which is particularly relevant during inflation. The quasi-exponential 
expansion stretches modes to vast superhorizon scales. Hence, a long wavelength approximation for the study of 
inhomogeneous spacetimes smoothed on scales larger than the horizon becomes particularly relevant. The resulting 
picture of the inhomogeneous universe is quite simple. Each point evolves like a separate homogeneous universe with 
slightly different values for the Hubble rate, scalefactor, scalar field etc. One can find nonlinear variables describing 
the inhomogeneity which do not depend on the choice of time hypersurfaces and which are exactly gauge-invariant in 
the perturbation theory sense. They are combinations of spatial gradients and are immediately connected with the 
linear gauge invariant variables 1Z and £ (see, for example, |5j and references therein). In this letter we would like to 
derive such variables describing the long wavelength curvature perturbation and find their evolution equations. The 
latter are essentially the same as those of linear perturbation theory. In the process we also elaborate on the separate 
universe picture which, so far, has only been used heuristically. 

II. THE LONG WAVELENGTH APPROXIMATION AND THE SEPARATE UNIVERSE PICTURE 

We start by assuming that we can smooth out all relevant quantities on scales much larger than the instantaneous 
comoving hubble radius ~ (aH)^ 1 . That is, for any quantity Q we will consider the quantity 



W is a window function rapidly decreasing for \x — x'\ > R and R is a smoothing scale larger than (aH)^ 1 . We 
assume that the Einstein equations act on the Q's. Of course this is not exact. Smoothing the full equations, i.e 
convolving the equations with the window function, is not obviously equivalent to the original equations acting on 
smoothed fields in the case of nonlinear operators. For example, if Q satisfies 



I. INTRODUCTION 



models. 




(1) 



D{Q) = 



(2) 



2 



then, in general, 

J d z x'D{Q{x'))W{\x - x'\/R) ± D(Q) (3) 

if D is a nonlinear operator. Here we assume that D(Q) = is an adequate description of the physics at long 
wavelengths. Then it is a reasonable approximation to say that we can drop those terms in the equations which 
contain more than one spatial gradient [2j. They are expected to be small compared with time derivatives. This will 
be the main approximation and will lead to a picture where different points of the universe evolve independently 1 . 

An interesting picture emerges when one also considers the quantum fluctuations of the various physical fields during 
an inflationary era. In such an era the comoving horizon decreases rapidly and the various modes that are exposed as 
they exit the horizon can be treated as classical stochastic fields |2^| . Assuming that the equations hold when acting 
on the Q's and writing the relevant Einstein equations as first order in time one sees that a time dependent window 
function results in extra terms in the equations. Indeed the time derivative of Q now has two terms, 

d -^= d ~i+j d 3 x'Q( X ')^W(\x - A/R) (4) 

The integration in (4) is restricted over modes close to horizon exit so Q(x') can be taken to be the quantum short 
wavlength field which has turned classical. Hence all evolution equations get augmented by a stochastic term. This 
is what has been used in the past to derive an effective langevin equation for the evolution of the inflaton field. In 
this letter we are concerned with the evolution of the spacetime after the inhomogeneity has been set up by quntum 
fluctuations. We develop a stochastic framework that takes into account gravitational perturbations and is gauge 
invariant in the long wavelength limit in [24(. 

Now, consider the ADM parametrisation of the metric 0, . One imagines spacetime to be foliated by spacelike 
hypersurfaces with a normal vector 

n = -N, n 2 = 0, n° = AT 1 , = N^N 1 . (5) 

Here, N is called the lapse function and N l the shift vector. The latter measures the deviation of the timelike lines 
which define the spatial coordinates from the integral curves of the normal to the time hypersurfaces. The metric 
takes the form 

.goo = -N 2 + 7«JV(JV ? , .go* = -N t , g V} = 7ij , (6) 

with the inverse 

g 00 = _ N -2 ^ g 0i = _ N -2 N i g ij = yj _ N ~2 N i N j ^ ( 7 ) 

where 7^ is the metric on the spatial hypersurfaces. The four functions N and N l parametrize the gauge freedom 
of general relativity and are arbitrary. The choice of gauge in the conventional language of perturbation theory is 
essentially the choice of N and N l for linearized perturbations. The way the spatial hypersurfaces are embedded in 
the 4D geometry is parametrised by the extrinsic curvature tensor 

Kij = ~[N i{j + N j[i + ^ i3 j. (8) 

Here a vertical bar denotes a covariant derivative w.r.t 7^. From now on we will use coordinate systems with shift 
vector N l — 0. For matter described by the energy momentum tensor T Ml/ we have the energy density 

£ = n^T^ = N- 2 T oa , (9) 

the momentum density 

Ji = -n^i = -N- 1 ^, (10) 



1 Particular solutions for long wavelength equations as well as an iteration scheme based on a g rad ient expansion has been considered by 
various authors in the past, see eg. [9|. A separate universe picture has also been discussed in lid in the context of perturbing solutions 
of the homogeneous equations w.r.t the constants that appear in them. Here, we are interested not in particular solutions but general 
properties of the nonlinear evolution of the long wavelength inhomogeneity codified by variables like (52). 
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and the stress tensor 

Sij — Tij. (11) 

A bar above a tensor will denote its traceless part 

Ki^Kij-^Kvi, K K\ -''A,,. (12) 

We can now write the Einstein equations for the above variables, dropping all terms which explicitely contain second 
order spatial derivatives and setting 7V*=0. The 00 and Oi Einstein equations are 

Kyle* - Ik 2 + ^£ = 0, (13) 
m P i 

R {\j - l K \i-^i = °- ( 14 ) 



mi 



The first will turn out to be equivalent to the Friedmann equation of the FRW cosmology and the second is usually 
refered to as the momentum constraint. The dynamical equations are 

dK% ( -_. 8tt -\ 

V - N \ KK '> - ^) • (16) 

where Sij = Sij — \S^ij with pressure S = S\ = -f tJ Sij and 1Z* is the Ricci tensor of the spatial hypersurfaces. Matter 
will obey the continuity equations T^ ; „ = which take the form 

d£ 1 
~dt 



— = NK(£+-S) + NK^Sij+N-\N 2 J%, (17) 

^ = NKJi (£Sf + Sj)N u NSi {j . (18) 

In the case of matter composed of several scalar fields <p A 1 as is relevant in scalar field driven inflation, the energy 
momentum tensor is 

= G AB d^ A d„<f> B - 9lxv (^G AB d x <f> A d x <l> B + V^j , (19) 



so that 



' G AB <j> A <j> B +V (</>), (20) 



2N 2 



Ji = ~G AB 4> A di4> B , (21) 

S V - Hi {^Gab^ B - V ) . ( 22 ) 

where we have dropped second order spatial gradients. So in the case of a collection of scalar fields in the long 
wavelength approximation (N 2 Ji)\j ~ and Sij ~ 0, up to second order spatial gradients. Actually, for a general 
perfect fluid 

= {£ + p)u^u u + pg^. (23) 

So, in the absence of pure vector perturbations on super-Hubble scales, itj = dif and the long wavelength approxi- 
mation gives S l j ~ pS l j, S l j ~ and J7i|j ~ 0. The set of relevant evolution equations then become 

dK 
dt 

dK) 



dt 
d£ 
~dt 

djj 
dt 





(24) 


NKK) , 


(25) 


nk(e + \s), 


(26) 


NKJi - (£51 + Sf)N u - NS>\. , 


(27) 
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which, along with the constraints 



K j 



£ = 0. 



16tt 
—Ji = 



form the basis of this approach. 

It will now be convenient to write the spatial metric as Q 

7ij = 7 1/3 7y 

where 7 is the determinant of the spatial metric, taking Det(jij) = 1. Now since Tr(7y) = 0, we have 



K, 



1 

'2N 



and 



K = 



_L1 



(28) 
(29) 

(30) 
(31) 
(32) 



The determinant 7(x, t) contains what in linear theory is usually called a scalar perturbation mode of the metric. 
7ij(x, i) contains another scalar, the vector and tensor perturbations. From (32) we see that if (x, t) can be interpreted 
as a locally-defined Hubble parameter since 7 1,/3 (x, t) = a(x, t) is a locally defined scalefactor. It can be explicitly 
shown that the system (24-29), a truncated form of Einstein's equations, forms a consistent set of equations, by which 
we mean that the evolution preserves the constraints. So an initial inhomogeneous configuration which respects the 
constraint (29), can be evolved with the evolution equations and the constraint always will be satisfied. Now observe 
that the only equation where spatial derivatives appear is the constraint (29) (and eqn (27) but this turns out to be 
irrelevant). Hence, after initial data which respect it have been specified, each point can be evolved individually. The 
relevant equations are essentially those of the homogeneous FRW cosmology, valid locally, modulo the K l j terms. 
This is exactly what has been termed in the past as the 'separate universe approach' (see, for example, |lj). ^ was 
explicitly stated in |2j in the context of general relativistic Hamilton-Jacobi theory. We see that such an approach is 
equivalent to the first order of a spatial gradient expansion. 

Although we can solve the equations numerically retaining the K l j terms, we can readily see that they are not 
expected to be important dynamically Q. Using (32) we can solve (25) to get 



K*j = C i i (x) 7 



-1/2 



(33) 



So, in the absence of sources, in an expanding universe (particularly for quasi-exponential expansion), the traceless 
part of the extrinsic curvature K % j decays rapidly. In most cases, therefore, it will be safe to ignore the K terms. 
Of course, in the presence of quantum fluctuations K is sourced mostly by the fluctuating 7 (see (31)). Ignoring K, 
equations (24-29) become 



dK 
dt 


l ^-N(G AB U A U B ) 
m P i 


— n A = 

dt 


NKU A - AT^ c n s H 


K 2 = 




diK = 


m P i 



NG 



. IB 



dV 



V 



(34) 
(35) 
(36) 
(37) 



where 



V = n ' 



(38) 



and r^ c is the Connection formed from Gab ■ Eqns (34-37) are exactly the same as those of a homogeneous cosmology 
apart from the momentum constraint (37) . Hence, the long wavelength universe looks like a collection of homogeneous 
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universes, each evolving according to the equations of FRW cosmology. The only information about inhomogeneity 
is contained in the momentum constraint (37) 0. 

In single field homogeneous cosmology, it is very convenient to parametrize the Hubble rate H in terms of the 
value of the scalar field leading to what is called the Hamilton Jacobi formulation. A similar thing can be done here. 
Observe that (37) is satisfied if K = K(cj) A ) and 

(39) 



yrA = ■■ L P l C AB 9K 

12tt d<t> B ' 



Substituting this into the Friedmann equation (36) gives the Hamilton Jacobi equation for a system of scalar fields 




R 2 _ ^ G AB 

12-k 



dK dK 24tt, 



-V. 



(40) 



"p1 



So, K is a function of the scalar fields only. Note that this is consistent with equation (34) which can now be written 

dK ■„ 



K = 



(41) 



Hence, given appropriate boundary conditions, a solution to equation (40) determines the state of the system com- 
pletely in terms of its position in the configuration space of the scalar fields. This result now holds for an inhomogeneous 
universe smoothed on supcrhorizon scales. 

Solving (40) analytically is not straightforward for a general potential. Explicit solutions for the case of an ex- 
ponential potential have been given in [2|. During inflation, it is usually a very good approximation to consider a 
slow-roll behaviour where the derivative terms can be neglected. Then (40) becomes the usual slow-roll relation 

if 2 = 2 -^V. (42) 

In this case the momentum constraint is equivalent to (36). Another case where an approximate solution can be found 
is the study of deviations of the scalar field about some stable point which are small compared to m p i 0- In any 
case, a solution to equation (40) describes the long wavelength system completely. An interesting point noted in 0, 
is the attractor property of solutions of (40). The latter contain a number of arbitrary parameters C A , reflecting the 
freedom to choose initial conditions of the field momenta (to specify the motion one must assign not only an initial 
value for </> but also at each point). Neglecting the K terms means that these parameters are spatially independent. 
Of course, when assigning initial conditions the momenta need not be strictly determined by the field values. So, in 
principle these constant parameters should be spatially dependent too. So suppose one considers a function K(cj) A , C), 
then 



2K™ = 


2 

"V qAB 


dc 


67T 






= 


m pl qAB 




12tt 












Ndt 



dK d dK 
d(b A d(b B ~dC 



dK d 



dK 



In 



dK 



dC 



dC 

where we have used (38) and (39). From (43) and (32) we see that 

dK 



dC 



oc a 



-3 



(43) 



(44) 



Thus the freedom to choose the momentum independently at each point, reflected in the freedom to choose constants 
in the solution of (40), becomes irrelevant very quickly 2 . 

From now on we will also be using the more familiar notation 

K{t,x) = -3H(t,x), 

7(*,x) = a 6 (i,x), (45) 
where a and H are the locally defined values of the scale factor and the Hubble rate. 



2 This is intimately connected to the fact that perturbation modes freeze when they exit the horizon 
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III. THE LONG WAVELENGTH CURVATURE PETURBATION 



In linear perturbation theory, a common variable used to characterize the perturbations is (5£ defined as 

6( = ^ - H 6 -i = $ - | (j, + H^j = j> - H 5 S (46) 

where ip is the scalar perturbation of the spatial part of the metric, V — 8a(t,x)/a(t). The quantity <5£ is gauge 
invariant and corresponds to the curvature perturbation in a time slicing that sets 5£ = 0. In the case of single 
field inflation, where there is one degree of freedom, on long wavelengths, 5( — 0. When there are many scalar fields 
present Q 

*C = ~7c~T — s^Pnadi (47) 
{£ +p) 

where Sp na d is the non-adiabatic pressure perturbation. In general, when there is a well defined equation of state 
p = p(£), 5p nac [=0- We will now derive a nonlinear generalisation of <5£ and eqn (47) for the nonlinear long wavelength 
case. 

Consider a generic time slicing and some initial time i/. We can integrate the local metric determinant up to a 
surface of constant energy density using (32) 

T 

Ino(t/,x) = lna(i,x) - / N(t' ,x)H(t' ,x)dt' , (48) 



ti 



where T(£(x)) is the value of t at the point x for £(x) = const. Taking the spatial derivative and using (32) once 
more we see that 

T 

ft[lna(t 7 ,x)] = U- ^d^lna^x) - J diiNit'^Hit'^dt' 

ti 

T 

= [di - ^dk\ In7(i,x) - J did t pna(f ,x)] dt' (49) 



which can be written 

dT 

9 s ;[lna(T,x)] = ft[lna(t,x)] - -^N(t,x)H{t,x), (50) 

or 

where we have defined Xi = di In a and T = T(£). Note that the gradient of the l.h.s of (50) is evaluated along the 
hypersurface £ — const. We see that with any choice of N, the particular combination of gradients appearing on the 
r.h.s of (51) always equals the gradient of In a on a surface where £ = const. So we are led to examine the variable 

NH 

Ci=Xi~ —d t £, (52) 

as a coordinate independent measure of the nonlinear curvature perturbation. Let us first see how this variable evolves 
with time. Taking the time derivative of (52) and using (26) we arrive at 

6 = ^ - 1 2 (mS - SdiS) , (53) 

or, by using (26) once more 
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This should be compared with (47). It is easy to see that if there is a well defined equation of state S — S(£), then 
the r.h.s of (53) is zero. Therefore, the long wavelength curvature perturbation, described by the vector Q, can evolve 
on superhorizon scales if and only if there exists a non-adiabatic pressure perturbation. This was proved in Q for the 
case of linear perturbation theory by making use only of energy conservation. The same applied here. Hence equation 
(28) shows that this result extends to the nonlinear case if we use the vector Q to describe the perturbation. 

The connection with the usual linearized gauge invariant variable £ is obvious from (52). Writing X — 111(7), 
£ = £{t) + 5£, 7 = 7(f) + <$7 and dropping terms quadratic in the perturbations, we see that 



<5C * --(^ ) (t,xo)5£ + C(t), ( r,r,) 



Q ^ di5(. 

with C(t) the background number of e-folds. 

The quantity Q was defined with respect to the energy density. Also the way we arrived at it was at best heuristic. 
We will now show that for any inhomogeneous spacetime scalar </>(i,x), the variable 

NH 

TZi =X % ^dit (56) 

<P 

is gauge invariant in the long wavelength approximation under mild conditions regarding the allowed coordinate 
transformations. Consider the transformation 

(t,x)-(T(t,x),X(i,x)). (57) 

Then, it can be shown that 2] 

Xi=f(x i ) + J jSf~ Q dT, (58) 

where / J is an arbitrary function independent of time. One can choose it for convenience to be SjX^ Q but one 
can also attach a physical meaning to such a choice as can be seen in figure 1. Related considerations in matching 
perturbed and homogeneous spacetimes are discussed in ^lj. We note that the transformation (58) holds between 
gauges with the shift Ni — 0. Assuming that the new time coordinate T is non-singular, that is, that the new time 
hypersurfaces do not get too 'wrinkly', we can discard the gradient of the integral as second order in spatial gradients. 
We then have 

Consider now a gauge with coordinates x^ where di<f> = and a transformation to a different gauge with coordinates 
X^. Then 

dX^dX* 
ox K ox 
~ -T k T t N^ + 6~ k k S l lTkI 

* ^JlTu, (60) 
7 (t,x)~7(T,X). (61) 



and so 



Taking the spatial derivative w.r.t x % we get 



dX 

For the scalar <j>(t, x) we have di<f> = 0, so in the new coordinates 



— In - —d 1 ~ — 9 ln~ 

dx % ^ dx l T ^ dx l dXi ^ 

d 

~ d l Td T lnj+ (62) 



OtTdrHT, X) + Qj^kT, X) = 0. (63) 
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FIG. 1: Two different slicings of an inhomogeneous spacetime based on two different time variables, T and T, have different 
normal curves denning the spatial coordinates. A particular patch of spacetime during inflation starts its life inside the horizon 
where it can be considered as homogeneous. Classical perturbations are generated when modes cross the horizon and freeze in. 
While subhorizon and homogeneous, there is a preferred time slicing in the patch to which all slicings and spatial coordinate 
choices should match. This fixes f 3 (x l ) = &\x % . This choice separates changes in the time slicing - which we want to consider 
as a gauge choice - from possible coordinate transformations of the homogeneous spacetime. 



a i ln 7 = ^-ln 7 -^— — 0. (64) 



Therefore 



Now, since the transformation we used was arbitrary (up to considerations for the smooth behaviour of the new time 
T), and noting that 

H = --±-d T l^ (65) 

and Xi = di In a, we conclude that the variable 

NH 

TZi =X t - -5-0i& (66) 
4> 

equals the gradient of the spatial metric determinant in a gauge where the scalar <j> is homogeneous; it is, therefore, 
a nonlinear gauge invariant measure of the perturbations. When the scalar is taken to be the energy density £ we 
recover Q of eqn (52). 

We would like to remark here that the use of spatial gradients to describe perturbations in a gauge invariant 
manner was first advocated in [l^. There, the use of the word 'gauge' refers to a choice of a correspondence between 
a fictitious background homogeneous spacetime and the real perturbed universe. So a gauge invariant quantity is one 
that does not change when this correspondence is altered, which is achieved when the quantity either vanishes or is 
a trivial tensor in the background |l3l Il4| . In this sense, any gradient is gauge invariant since it must vanish on the 
homogeneous background. Here, however, by a choice of gauge we mean a choice of the lapse function N and the 
shift vector iV, |5j, so we look for quantities such as (66) that are invariant under a change of time slicing as well. 

Examples 

When the scalar in (56) is taken to be the inflaton it can be easily verifield that TZi is also conserved. Using eqns 
(34-37) one can show that 

Ki = . (67) 

When there is more than one scalar field present, however, the situation is different. It is well known in linear 
perturbation theory 0, ITU Il6| that perturbations in multi-scalar field models can source the curvature perturbation 
<5£ on long wavelengths, making the time evolution of the latter possible. Let us see how this effect is very easily seen 
in the formalism developed in the previous section without the need to resort to linear perturbation theory. Consider 
hypersurfaces with di£ — 0. From (28) with the K terms omitted we see that 

£ ti = H,i = Ji = G AB Tl A d^ B = 0. (68) 

The fact that = also means that 

S ii = 2G AB U A d i U B . (69) 
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So, if there is a single field, 

£ A = => H ti = => d l( j) = => S,i = 0, (70) 
since 9^11 = 0. Hence Q remains constant. On the other hand, for many fields, we have 

e ti = o => = o => G AB n A a^ B = o. (71) 

So, generally in this case, we may have di<j) A ^ 0. Therefore, from (69) diS ^ in general and Q can evolve on 
superhorizon scales, Q ^ 0. 

In general, for a system of scalar fields governed by the Lagrangian 

C = \G AB d^ A d^ B - V, (72) 



the variation of Q from (54) is given by 



Q = V A 



1 N d,6 A - s A d t ' 1 



(73) 



dt\G BC U B U c J lT r l \G BC U B U c / 

In the multi field case it is convenient to consider the variables 

Q A = \l A X t - Hdi4> A , (74) 

which are also invariant under changes of time hypersurfaces in the long wavelength limit. They are not conserved 
even in the case of a single field. In particular 

Qf = tl A K A +T1 A K A , (75) 

where no summation is implied. It should be noted here that it is possible to define conserved variables for each field 
if and only if the potential is of the form ^ V (4> A ) or if one field dominates completely over the other. In terms of 
the Q A, s the curvature perturbation is given by 

= G AB U A Q B 

^~G CD Tl c Tl D ' 1 ' 

Consider now the following simple two-field example which can be considered a toy model for hybrid inflation. Initially 
the first field a is stuck slightly displaced from the bottom of a valley while the other field <j> evolves and drives inflation. 
In this case, we have essentially a single field model and £i is conserved. With W ~ 0, we have 

„(/) „ H ( „ , . Gd 



When the <f> field reaches a critical value, the trajectory turns abruptly in the direction of variations in a which can 
now dominate until inflation ends. For the period of rolling in the er-direction we have an extra contribution to Q 

+ (78) 

Note that for a non-diagonal metric G AB the isocurvature cr-field can contribute to the curvature perturbation even 
for a straight trajectory m contrast to what happens when G AB = S AB Also, if during the rolling along 

the a direction we have a second period of inflation, the a field can contribute to the curvature perturbation since 

In the case with G AB = 6 AB , we can easily find the following evolution equations for the IZi's in the case of two fields 

Uf = -^-ATI ff [dio - ) (80) 

m P i V 11 / 

and 

4tt „ n, 
5 



Tlf = — iVn di<t> - -?-di<T ) . (81) 

m r,l V W<T 
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For many general fluids one can write the curvature perturbation as a weighted sum over components |25j 

d = (82) 

If the components are non interacting, each will obey an evolution equation like (26), so each is separately 
conserved. Hence one can see that the time evolution of the total curvature perturbation is determined by the ratio 
of the energy densities of the various components. For example, starting with C,f — we get 

Ci - T^tr^? (83) 

OA + OB 

This evolution of the long wavelength curvature perturbation in multifield models is one of the possible mechanisms 
that could potentially produce significant non gaussianity. Perturbations from the a direction can be significantly 
non-gaussian and now there is the possibility that these fluctuations can contribute to the curvature perturbation 
at later times [19J, |2fJ . This can be seen for example from equations (80) and (82). Other mechanisms relating 
non-gaussianity with a second scalar field have been proposed |2lj |. We will return on the issue with more details in 

One can also easily see that there is a particular form of a multi-scalar field potential which will conserve With 
= we have G ab^1 a di<f> B = 0, so to get — we need (see (69) and (39)) 

Gab&IW* = ^GabIIM (g bc ^\ = . (84) 



Thus, if 

GAU lTIr = —<t> A => iTA = —Gab^ b , (85) 



., AC dH _ A A _^ dH _ _>^ r u 
m p i d(j) A m pl 



with A some dimensionless number, then S t i = 0. From (40) 



H2 = 'ik GAB w*w + ^ Vi '- 



8tt 
5 



So, with 

A 1 



H = —^-Gab^^ + /im pl , (87) 
m p i 2 



we get 



V{4>) = -Im 2 4> A 4> a + ^ A ^A + qm 2 p ^ 



pl 

with q another dimensionless parameter and M 2 can have either sign. So, for a general potential containing both 
quartic and quadratic terms with the various constants related as in (76), the long wavelength curvature perturbation 
does not evolve in time. One should keep in mind that the above results only apply when anisotropic stresses are 
negligible on superhorizon scales and therefore one can safely neglect the terms. 

We emphasise again that for the derivation of all of the above no linear approximation has been made. The 
equations are exact and hold at every point of the inhomogeneous spacetime in the long wavelength limit. For more 
uses of this formalism we refer the reader to |24|. 



IV. SUMMARY 



We have shown that particular combinations of spatial gradients are good variables to describe the curvature 
perturbation in a long wavelength limit. We also elaborated on the intuitive picture of the 'separate universe approach' 
where, given initial data which respect the energy and momentum constraints, one can view the long wavelength 
universe as a collection of universes each evolving independently. Such a picture is a simple way to understand 
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nonlinear evolution on large scales and, in particular, to determine the conditions under which the long wavelength 
curvature perturbation can vary. We believe this is a remarkably straightforward and physically intuitive approach 
which should find much wider application. 

Note added: While this paper was undergoing final revision, independent work by Lyth and Wands appeared on the 
archive |25|. This addresses some related issues on ^-conservation using linear perturbation methods. 
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